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A  BOUNDARY  SCATTERING  STRENGTH  EXTRACTION  ALGORITHM  FOR 
THE  ANALYSIS  OF  LONG-RANGE  REVERBERATION  DATA 


1.  INTRODUCTION 


This  report  describes  the  Scattering  Strength  Extraction  Algorithm  (SSEA)  model,  devel¬ 
oped  to  extract  boundary  acoustic  scattering  strengths  from  measured  underwater  acoustic  re¬ 
verberation  time  series.  Ocean  surface  scattering  strengths  have  been  deduced  from  short-range 
data  with  the  results  of  Chapman  and  Harris  [1]  perhaps  being  the  most  well  known  and  widely 
used.  The  geometry  of  short-range  experiments  is  simple  enough  that  scattering  strengths  may 
be  extracted  by  combining  the  results  of  ray  traces,  beam  patterns,  and  geometric  calculations 
in  a  sonar  equation  analysis.  For  long-range  data,  the  surface  is  not  uniformly  ensonified  and 
there  are  multiple  ray  paths.  These  two  complications  make  it  impractical  to  extract  scattering 
strengths  from  long-range  data  without  the  use  of  a  comprehensive  computer  model.  The  SSEA 
model  allows  us  to  extract  scattering  strengths  from  measured  long-range  reverberation  time  series 
and  will  enable  us  to  corroborate  the  scattering  strengths  deduced  from  short-range  measurements 
with  those  from  long-range  data.  It  is  desirable  to  corroborate  with  long-range  data  for  several 
reasons: 


•  to  demonstrate  that  the  short-range  data,  typically  collected  by  using  impulsive  broadband 
signals,  can  be  used  to  predict  reverberation  from  all  signal  types; 

•  to  extend  the  available  scattering  strength  data  base  with  regard  to  enviromental  conditions 
and  signal  types. 

2.  REVERBERATION  EQUATIONS 


2.1  Short-Range  Reverberation 


by 


As  shown  in  Fig.  1,  the  calculation  of  reverberation  levels  for  short-range  geometries  is  given 


R(t)  =  ^2Q(t-Tou t„  -  Iback.  )6s  W’s. ,  <t>sn IMV’r,  ,  4>iu ) 


S (0out„  i  ^back,  i  l/O  An 


Lout„  Lbackn 


]■ 


(1) 


Here  Q(t)  is  the  source  envelope,  and  the  summation  over  n  represents  an  integration  of  scattering 
contributions  from  each  surface  area  element.  In  the  case  of  long-range  reverberation  it  js  necessary 
to  add  an  additional  sum  over  ray  pairs.  In  this  case,  however,  there  is  only  one  incident  and  one 
scattered  eigenray  to  each  surface  patch.  The  angles  ips*  and  <ps„  are  the  horizontal  and  vertical 
launch  angles  at  the  source,  and  65  is  the  relative  beam  power  in  that  direction.  The  same  relation 
holds  for  the  receiver  ray  arrival  angles  ipR*  and  and  for  the  relative  receiver  beam  power 
6r.  The  horizontal  bistatic  angle  i/>  is  the  difference  between  the  azimuthal  angles  of  the  incident 

Manuscript  approved  July  18,  1989. 


1 


PITRE,  FROMM,  AND  F RANCHI 


Fig.  1  -  Short-range  geometry 


and  scattered  rays.  We  have  included  the  horizontal  bistatic  argument  in  the  scattering  strength 
function  for  completeness  in  our  discussion  although  the  cases  that  we  ultimately  consider  are 
independent  of  V’-  For  each  n  in  the  surface  integration,  the  nth  surface  patch  of  area  An  is 
uniformly  ensonified  and  contributes  reverberation  characterized  by  specific  values  for  the  travel 
times  t;  the  transmission  losses  l/Lout  and  1/1/ back!  and  the  grazing  0outn  and  backscatter  angles 
^back*  •  Geometrical  spreading  from  source  to  surface  is  represented  by  l/Lout,  and  spreading  from 
surface  to  receiver  is  represented  by  1/Xback-  We  interpret  R(t)  to  be  the  multi-ping  average  of 
the  pressure  squared,  where  t  represents  time  with  respect  to  ping  transmission. 

Experiments  to  determine  backscattering  strengths  are  usually  designed  so  that  the  source 
and  receiver  beam  patterns  and  a  finite  ping  duration  selects  a  small,  uniformly  ensonified  surface 
patch  for  which  0out  =  #back-  In  this  situation,  it  is  only  necessary  to  consider  one  term  in  the  sum 
over  areas.  By  taking  10  log10  of  the  summand  (contribution  from  a  single  surface  area  patch)  of 


Eq.  (1),  we  get  the  reverberation  sonar  equation 

RL  =  SL  —  AG  +  SS  4- 10  logjo(i4n)  +  TL0 ut  +  TLback  >  (2) 

where 

SL(t)  =  10  log10  [Q(t  —  Tout  —  Tback)]  (3) 

AG  =  10  iog10  [MV’S.  <£«)]  (4) 

SS  =  10  loglo[S(0)]  (5) 

TTout  =  10  logio  [1/Tout]  (6) 

TLback  =  10  l°8l0  (1/ I'back]  •  (7) 


2 


NRL  REPORT  9142 


2.2  Long-Range  Reverberation 

The  calculation  of  long-range  reverberation  levels  is  similar  to  that  of  short  range,  but  with 
the  addition  of  a  summation  over  contributions  from  source-receiver  multipath  scattering  pairs. 
Figure  2  illustrates  three  possible  ray  paths  out  to  the  nth  surface  patch  in  the  first  convergence 
zone.  The  multipath  summation  J  for  these  paths  produces  nine  contributions  to  the  reverber¬ 
ation.  In  general 


m  =EEW- T".«  -  TnjM+S.J+S.JbR&R.sAR.j) 


S{0n,iy^nji‘4>ij)  -An 

Ln,» 


(8) 


The  index  i  denotes  a  source  ray,  and  the  index  j  denotes  a  receiver  ray.  Equation  (8)  is 
the  basis  of  an  extensive  long-range  reverberation  model,  the  Range-dependent  Active  System 
Performance  Prediction  Model  (RASP),  which  generates  reverberation  time  series  that  are  used 
in  active  system  performance  prediction.  Reference  2  describes  this  model  and  its  application. 

3.  SCATTERING  STRENGTH  EXTRACTION  APPROACHES 

An  initial  approach  to  extracting  scattering  strengths  by  use  of  Eq.  (8)  employs  the  fact 
that,  at  surface  convergence  zones,  the  energetic  acoustic  rays  usually  span  a  small  angular  range. 
Thus,  if  one  assumes  that  at  t  an  average  grazing  and  backscatter  angle  can  be  used,  the  resulting 
equation,  in  dB,  is 

RL(t)  =  SL  +  SS(0ave)  +  10  log}Q  Roo(t)  ,  (9) 

where  Roo(t)  is  a  model  prediction  with  a  0  dB  source  level  and  a  0  dB  scattering  strength. 
Comparing  RL(t)  with  the  corresponding  experimental  time  series  yields  average  values  for  the 
scattering  strength  for  the  nominal  grazing  and  backscatter  angles.  This  process  was  used  with 
some  degree  of  success  on  both  artificial  and  simulated  data.  However,  surface  scattering  strengths 
are  very  sensitive  to  grazing  angle,  and  a  more  refined  extraction  process  is  needed  to  retain  the 
grazing  and  backscatter  angle  dependence  in  the  results. 
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The  first  refinement  considered  was  an  inversion  of  the  matrix  equation  that  results  from  the 
discretization  of  t  in  Eq.  (8).  This  method  was  not  attempted  for  several  reasons  relating  to  the 
number  of  time  samples  needed  to  generate  an  invertible  equation. 

Ultimately,  it  was  found  that  an  extension  of  the  original  method  to  a  two-parameter  fit  (as 
opposed  to  a  single  parameter  0  dB  fit)  was  tractable  as  well  as  sufficiently  general  to  be  useful. 
A  parametric  form  for  the  scattering  strength  function  is  selected,  and  its  parameters  are  varied 
to  obtain  an  optimum  fit  to  the  measured  reverberation  time  series.  The  critical  aspects  of  this 
method  are  the  selection  of  a  physically  reasonable  scattering  strength  function  with  a  minimum 
of  fitting  parameters  and  use  of  the  optimization  algorithm. 

4.  SELECTING  A  FUNCTIONAL  FORM  FOR  SCATTERING  STRENGTHS 

We  have  chosen  the  following  two-parameter  scattering  strength  function: 

5(^backt  ^out)  =  a[sin(0out)  sin(0back)]7-  (10) 

The  scattering  strength  parameters  a  and  7  are  implicit  functions  of  sea  state  and  source  fre¬ 
quency.  Parametric  forms  for  a  and  7  as  functions  of  sea  state  and  frequency  can  be  determined 
by  applying  SSEA  to  data  taken  at  many  different  frequencies  under  a  variety  of  sea  state  condi¬ 
tions.  There  is  no  horizontal  bistatic  angle  dependence  to  our  scattering  strength  function  at  this 
point.  Applications  of  the  SSEA  method  to  horizontally  bistatic  source  receiver  configurations 
will  require  the  use  of  parametrized  scattering  strength  functions  with  a  horizontal  bistatic  angle 
dependence.  The  methods  and  formulations  that  are  introduced  below  change  very  little  in  their 
details  if  scattering  strength  functions  with  more  parameters  are  introduced;  however,  the  level  of 
computational  effort  can  increase  dramatically. 

Our  chosen  form  for  the  scattering  strength  function  is  consistent  with  the  small-angle  form  of 
the  Raleigh  perturbative  scattering  strength  for  an  ocean  surface  described  by  a  Pierson-Moskowitz 
spectrum: 

8 1  x  10~3 

*5 Rayleigh  (^outi  $bac  k)  ~  (^out^back)  •  (U) 

It  is  also  consistent  with  the  Chapman  and  Harris  empirical  fit  to  surface  backscatter  [3]  derived 
from  monostatic  (0Out  ~  #back)  data 

a(f,v)  (0P(/’O)  ,  (12) 

where  /  and  v  are  the  acoustic  frequency  and  wind  speed,  respectively.  Thus  the  choice  of  a 
two-parameter  fit  of  the  form  of  Eq.  (10)  for  the  scattering  strength  function  has  its  basis  in  both 
theory  and  experiment. 

A  three-parameter  form  of  the  scattering  strength  function,  <*(0Out#back)7  +  was  examined 
because  of  its  possible  relation  to  scattering  that  originated  from  a  subsurface  layer  of  volumetric 
inhomogeneities.  The  resulting  data  fits  were  not  acceptable  at  low  angles  since  negative  linear 
scattering  strengths  were  predicted. 

Linear  forms  for  the  scattering  strength  function  like  ad  4-  7  were  also  considered  in  the 
expectation  that  the  scattering  strength  function  would  be  well  enough  approximated  when  the 
scattering  angles  contributing  the  the  reverberation  vary  over  a  small  range.  Moreover,  a  linear 
scattering  strength  form  simplifies  the  least-squares  inversion  process  described  in  the  next  section 
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so  that  it  reduces  to  a  matrix  inversion.  The  angle  domain  over  which  linearly  extrapolated 
scattering  strengths  are  useful  was  found  to  be  too  narrow  for  practical  applications.  Rather, 
the  logarithm  of  the  scattering  strength  behaves  linearly  over  a  large  range  of  the  logarithm  of 
scattering  angles. 

5.  SCATTERING  STRENGTH  EXTRACTION  ALGORITHM 


To  determine  the  scattering  strength  parameters  a  and  7,  we  use  a  least-squares  approach. 
The  RASP  reverberation  model  is  used  to  generate  the  average  reverberation  time  series  R(t,  a,  7) 
resulting  from  a  given  choice  of  the  scattering  parameters.  The  error  that  we  mimimize  is  the  time 
integral  of  the  squared  difference  between  this  trial  series  and  the  actual  data  R^it): 


E(a,  7) 


/  1  \  rtj  [R(t,a,7)  -  RAf(t)]2 

\t2-tJJu  Rlf(t) 


(13) 


It  is  important  to  note  that  the  integrand  is  actually  the  square  of  the  relative  error.  An  abso¬ 
lute  error  condition  was  used  in  the  initial  implementations  of  the  method.  Because  scattering 
strengths  and  transmission  losses  range  over  several  orders  of  magnitude  in  the  integration  window, 
it  is  not  always  possible  to  determine  the  scattering  strength  parameters  with  an  absolute  error 
minimization. 


The  error  is  minimized  by  choosing  a  and  7  so  that 


dE(a,  7) 

da 


and 


dE(a,  7) 

07 


Because  the  scattering  strength  is  linear  in  a,  and  the  reverberation  is  linear  in  the  scattering 
strength,  we  can  solve  the  extremum  condition  for  a, 


and  then  we  are  left  with  the  nonlinear  equation  for  7: 


(14) 


q(7)R(t,l,7)  ~  RmW 

RlfW 


0R(t,  1,7) 

07 


dt  =  0. 


(15) 


Note  that  R(t,  1, 7)  and  dR(t,  1, 7)/07  are  equivalent  to  model  reverberation  time  series  predictions 
with  the  scattering  strengths  given  by  (0Out,0back)7  and  ( 0OU  t  ^back ) 7  In  ( #ou  t  #back ) ,  respectively. 
There  may  be  many  solutions  to  Eq.  (15)  in  a  long-range  multipath  case,  but  usually  one  solution 
is  well  isolated  from  the  rest,  which  is  commensurate  with  the  expected  range  of  values  of  scattering 
strength  for  the  appropriate  wind  speed  and  frequency. 


6.  RESULTS  FROM  ANALYZING  THE  SIMULATED  DATA 


Figure  3  shows  a  model  prediction  of  surface  and  bottom  reverberation  for  a  typical  deep 
ocean  enviroment.  The  signal  was  a  250  Hz  CW  with  a  pulse  duration  of  10  s,  and  Chapman  and 
Harris  surface  scattering  strengths  for  a  17  knot  wind  were  used.  Note  that  outside  of  the  surface 
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convergence  zones  the  reverberation  is  dominated  by  bottom  backscatter  so  that  our  analysis  of 
surface  backscattering  strengths  must  be  confined  to  regions  near  the  convergence  zone.  Figures  4 
and  5  show  the  results  of  applying  SSEA  to  several  small  time  windows  on  and  about  the  second 
convergence  zone  peak  (at  150  s).  The  behavior  of  the  time  senes  is  quite  different  in  each  of 
these  windows.  Reverberation  in  Fig.  3  is  roughly  montonically  increasing  in  the  earliest  window 
(140  to  150),  montonically  decreasing  in  the  latest  window  (150  to  160),  and  rising  and  then 
falling  behavior  in  the  largest  window  (140  to  160).  The  near  constancy  of  the  scattering  strength 
parameters  extracted  from  these  three  windows  indicates  a  degree  of  robustness  in  the  ability  of 
SSEA  to  deconvolve  the  complicated  convergence  zone  transmission  loss  behavior  and  the  simple 
monotonic  scattering  strength  function  from  the  reverberation  time  series. 

7.  NOTES  ON  THE  APPLICATION  OF  SSEA  TO  SURFACE  REVERBERATION 

The  example  application  given  above  demonstrates  the  efficacy  of  SSEA  under  ideal  condi¬ 
tions.  In  this  section  we  indicate  a  few  of  the  outstanding  issues  that  have  arisen  in  its  application 
to  actual  data. 

Because  the  convergence  zone  time  intervals  in  which  the  surface  reverberation  dominates 
over  bottom  reverberation  are  short,  it  is  important  that  time  averaging  in  the  data  preprocessing 
be  kept  to  a  minimum.  Preprocessing  time  averaging  windows  that  are  a  significant  fraction  of  tlie 
convergence  zone  width  have  Wen  found  to  significantly  alter  the  scattering  strength  parameters, 
even  when  the  averaging  process  was  emulated  by  the  reverberation  calculation. 

Depending  on  the  signal  type  it  has  in  some  cases  been  enlightening  to  filter  the  reverberation 
time  series  into  Doppler  bins  and  compute  scattering  parameters  for  each  of  these  bins  separately. 

It  is  important  to  note  that  SSEA  assumes  that  any  differences  between  the  experimental 
data  and  the  model  prediction  are  due  to  the  scattering  strength.  If  a  significant  amount  of  time 
has  elapsed  between  the  pings  used  to  produce  a  ping-averaged  set  for  analysis,  then  variations 
in  parameters  such  as  the  bottom  loss,  the  sound  speed  profile,  or  the  source/receiver  beam 
pattern  orientations,  may  not  be  accurately  represented  by  the  modeling.  These  inaccuracies  in 
the  modeling  will  manifest  themselves  in  the  extracted  scattering  strengths.  For  this  reason  it  is 
strongly  recommended  that  pings  be  back  to  back  and  in  nonoverlapping  frequency  bands  to  get 
useful  ping  averaging. 

8.  SUMMARY 

We  have  developed  a  long-range  scattering  strength  extraction  algorithm.  It  has  been  suc¬ 
cessfully  tested  on  modeled  reverberation  time  series  generated  with  specific  scattering  strength 
functions  and  has  been  shown  to  recover  the  original  scattering  function. 

The  SSEA  model  represents  a  new  capability  for  extracting  detailed  information  on  the 
scattering  processes  from  the  analysis  of  long-range  reverberation  data.  It  is  still  being  refined  to 
provide  improved  accuracy  and  to  account  for  more  details  of  the  signal-dependent  signal  processing 
methods  used.  A  short-range  bistatic  model  is  also  under  development.  This  report  focuses  on 
the  use  of  SSEA  to  extract  surface  backscatter  parameters,  but  it  is  equally  applicable  to  the 
extraction  of  bottom  backscatter  parameters.  In  fact  Lambert’s  Law  for  bottom  backscatter  has 
the  same  form  as  our  chosen  parametric  scattering  strength  function. 


6 


BEAM  POWER  IdB  re  1  M P®2> 


NRL  REPORT  9142 


100 

90 

60 

70 

60 

60 

40 

0  50  100  150  200  250  300  350  400 

TIME  (S) 

Fig.  3  -  Reverberation  using  a  Chapman  and  Harris  surface  backscatter  function  for  250  Hz  and  17  knots 
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Fig.  4  -  Comparing  Chapman  and  Harris  scattering  strengths  with  those  extracted  by  SSEA  for  various 
time  intervals  that  bracket  the  peak  of  the  convergence  zon*' 
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Fig.  5  -  Comparing  Chapman  and  Harris  scattering  strengths  with  those  extracted  by  SSEA  for  various 
time  intervals  preceding  and  following  the  peak  of  the  convergence  zone 
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